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Abstract 



The exact reduced dynamics for the independent oscillator model in the RWA approximation at 
zero and finite temperatures is derived. It is shown that the information about the interaction and the 
environment is encapsulated into three time dependent coefficients of the master equation, one of which 
vanishes in the zero temperature case. In currently used optical cavities all the information about the 
field dynamics is contained into two (or three) experimentally accesible and physically meaningful real 
[ functions of time. From the phenomenological point of view it suffices then to carefully measure two 

' (three) adequate observables in order to map the evolution of any initial condition, as shown with several 

L~:> ' examples: (generalized) coherent states, Fock states, Schrodinger cat states, and squeezed states. 



1 Introduction 



' Measuring the time development of the entanglement process of a system coupled to its environment is a most 

^ I remarkable achievement and a challenging goal. The reason for this is that the entanglement process is a 

unique and typical quantum feature. Several attemps, both on the theoretical as well as on the experimental 
side have been recently made|l|, |^. In the particular case of high Q optical cavities a direct measure of 
the decoherence process has been given and suggestions of experiments with essentially the same set up 
^ have been made on how to directly measure the Wigner function of the initially correlated field produced in 

■ the cavity||, ^. We show that all necessary information to construct the time development of any Wigner 

function (or any system's density operator) can be obtained by a precise measure of three quantities as 
function of time: the average photon number and two orthogonal field quadratures. 

Recently the exact master equation for quantum Brownian motion in a general environment has been 
derived using both path integral techniques Q and the tracing of the evolution equation for the Wigner 
function||, |^. We closely follow the later approach to derive the exact master equation for the oscillator 
independent model in the RWA at zero temperature and with a factorized initial condition. The hamiltonian 
of the model is 

H = huj{a)a + 1/2) + Ti'^ujk{alak + 1/2) + Ti'^Ck{a)ak + a\a). (1) 

k k 

Here we present the solution of the initial value problem, which have been solved in the Heisenberg picture 
in Refs. p, p|. The solution allows for an easy visualization in contrast to the model without the RWA 
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approximation. Besides its intrinsic interest as an exactly soluble model, the hamiltonian may be useful 
in treating leaking Bose-Einstein condensates ]lO| , in materials with modified dispersion relation s or in 
any case of non-ohmic strength function, where the Born-Markov approximation is not adequate]lO|. 
We assume that a.t t — the total density operator is given by 

where the subscript S for system refers to the main oscillator. The bath, i.e. the set of oscillators labelled 
by fc, is initially in thermal equilibrium at inverse temperature p. At zero temperature, the tensor product 
of the vacuum of the main oscillator with the vacuum of the set of oscillators is the ground state of (^. 

2 The Exact Master Equation 

It is well known that for quadratic hamiltonians the Wigner function satisfies the classical Liouville equation. 
To obtain the classical hamiltonian corresponding to , one uses the correspondence rule , where 



and discards the zero energy contributions. We use /x = 0, 1, 2, .., k — 1,2, .. and oq = a, luq = lo. In general 
greek subindices denote non negative integers while latin sub indices denote positive integers. Using these 
conventions, the classical Liouville equation 

dW{a^,a*^,t) _ 1 dH dW 1 dH dW 
dt ih ^ doLu da* ifi ^ da* dan ' 

where W^(a^,Q;*,f) is the Wigner function in the quantum case and the probability density function in the 
classical case, can be written as 

dW{af„al,t) dW dW ^ ,dW ^ dW 

k " k 

k k k k 

The initial condition Eq. (|^) in the language of Wigner functions is 

W{a^,a*^,t^ 0) = W\a^,a*^) = W^s{o^,a*)W%{ak,al) = TyO(a, a*) [] ^^fee-2t-M?^-'=«"'=«^ (6) 

fe 

where the Nk are normalization constants. Integrating Eq.(|^) over the bath variables we get 

dW(a,a*,t) dW dW dG* dG 

r — —lua — h lua— I— h (7) 

dt da* da da* da 

with 



and 



W{a,a*,t) = J i^dakdalj W{a^,a*^,t), (8) 
G{a,a*,t)^ J ^Y[dakda*^'^CkakW{a^,a*^). (9) 
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As in [Q it is easy to show that G{a,a*,t) can be written in terms of W. We notice that for quadratic 
hamiltonians VF(a^,Q!*,i) = VF°(a^(— /:), a* (— i)), where is the solution of the classical equations of 

motion. If we define d{t) = {ao{t), ai{t), . . .) and denote its transpose by d^{t), we have 

d^{t)^U^A{t)Ud^{0), d*{t)^{d{t))*, (10) 

where U and A are unitary, and A is diagonal. Taking the Fourier transform of G, and changing variables 
from {q!^(— t), a* (— t)} to {a^t(O), a* (0)}, with unit Jacobian, we obtain 

G(K,/t') = /"[]da^(0)da;(0)e*"S.'''^(*)"-(°)e'"'S.P^(*^""(°) 

xY,qAt)aA0)Wlia^*\0))WUa^*\0)\ (11) 

with {pij,p1, Qij} time dependent parameters. From Eq.(pl|) it is easy to see that the multiplication by ao = a 
is equivalent to a derivation with respect to fc, plus terms corresponding to multiplication by ak, up to time 
dependent coefficients. These last terms, as can be seen from (^, correspond to derivations w.r.t. a^, which 
in turn, are equivalent to multiplication by fe', as shows a simple integration by parts. Taking the inverse 
Fourier transform we^btain a multiplication by a and a derivation w.r.t. a*. Thus, observing that the 
Fourier transform of VF(q!, a*), W{k, k!) is given by 

W{k,h') = / []dQ^(0)da;(0)e*"S„P''(*)"-^We^-'E„P:(*)<(o) 
xW^°(aW(0))W^°(aW(0)), 



we obtain 



with time dependent functions Y, Z . Therefore the Wigner equation can be written as 

dW{a,a\t) .. ,J ,dW dw\ ^~ 

^ ' = -i{uj + 5)\a*—^a—\+2\W 

at \ da* oa I 

+A fa*— + + A'-^^ (13) 

I da* da f dada* ' 

where we have set iY = X + iS and iZ — iZ* = A'. All of A, 6 and A' are real functions. By comparing the 
system of equations found from both (I) and dH) we get 

{X + iS){a) = i^Cfc(afe), (14) 

k 

{X + iS){a^) = i^Ckiaak), (15) 

k 

-2X'{aa*) = i^Ck{aal) - i'^Ck{a*ak). (16) 

k k 

We know that the solution of the Heisenberg equations can be written as follows 

a{t) = r^{t)a{0) +J2lk{t)ak{0), (17) 

A: 

ak{t) - Vk{t)a(t)+Y,lkiai{0). (18) 



3 



Using the above solution of the Heisenberg equations, and the fact that aU of the first and second (symmetric) 
moments involving bath operators are zero, with the exception of {{al,,ak}) = 2nk + 1 = 2 coth{'hu!k(3/2), 
we obtain 

X{t)+tS(t) - iY,CkVk{t) (19) 

k 

x'{t) = J2(^k{nkm + l){tjn*ki-ntiH). (20) 

kl 

For reasons that will be clear soon, we write the diffusion coefficient X'{t) as X{t)+e{t, f3), with lim^^oo ^U, (3) — 
0, as shown in section ^ At zero temperature, since the tensor product of vacua is the ground state of (|l|), the 
corresponding (reduced) Wigner function PFs(a, a*) should be a stationary solution of the Wigner equation. 
When this condition is applied to Eq.(|T^), we obtain A = A'. If they were not equal it would imply the non 
existence of an exact master equation, as seems to be claimed in Ref. [p^ . However, this is not the case, as 
we show next. 

It is not hard to show that the operator equation for the system's reduced density operator equivalent to 
the Wigner equation (RSh is 



= — Ihiuj + S)a^a, p] + (A + e)(2a • — a^a • — • a^a)/9 + e(2a^ •a — aa^ • — • aa^)p = C(t)p(t), (21) 
dt in 

where the usual dot superoperator convention has been used. The usual Born-Markov RWA master equation 
is of this form with constant coefficients |[T^. Some results can be obtained at once from (|T]): premultiplying 
by a and taking the trace we get 

.(-.(. + . )-A)«, (22) 

which can be immediately solved to give 

ait) = exp(-il7(i) - A(t))a(0), (23) 

with ^ ^ 

n{t)^ I dT(w + (5)(T), A(i) = / drXir). (24) 



Note that this result is independent of e, i.e., it does not depend on the temperature. Premultiplying ( pl| ) 
by a^a, and taking the trace we get the following differential equation 

— (at a) (i) = -2X{a)a) (t) + 2t (25) 
dt 

with the solution 

(ata)(i) =exp(-2A(<)) (at a) (0)+A/'(<) = (at a) (t;/3->oo) + 2 exp(-2A(t)) \ dre(T) exp(2A(T)), (26) 

where it is evident that M[t) vanishes in the zero temperature limit. Contrary to the exact equations found in 
1^, 1^, 1^ We thus have the following interpretations for the real functions that appear in the master equation: 
(5(t) is the instantaneous frequency shift, A(t) is the instantaneous energy rate of change at zero temperature 
and tit) is the instantaneous energy rate of change at finite temperature but with the system in the vacuum 
state. Moreover M{t), which is related to both e{t) and S{t) is the mean number of excitations when the 
initial state was the ground state. 
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3 The evolution superoperator and some initial states 

We can use Lie algebraic niethods[0] to find the evolution superoperator U. Indeed, we can verify that the 
superoperators M = a^au, V ~ •a'a, J ^ a* and = • a form an algebra, 

[x,7']-o, {M.,j\ = - J ^{v,j\, [^^,7^] = -7^ = [7',7^]. (27) 

Thus, we can assume that U (t) = we'^^e^^e^^'e^"^. Deriving this expression, using the formula ex\){xA)B exp(— a;A) 



B + x[A, B] +x [B, [_B, A]]/2! + ... and the commutation relations (27), comparing coefficients in the equation 



dUldt = C{t)U{t), and solving the resulting differential equations we obtain 

exp(-2A(t)) 

'^'^ = ^- i+m ■ 

Let us suppose that p satisfies the equation dp/ dt = C{Xi*^ 'Xf, t)p{t), where the Xi are operators (notice 
that £ is a general linear superoperator), and that p can be written as Up'U^^. Then, p' satisfy the equation 
dp' /dt = C'p'{t). If, moreover, we choose U = exp((Ta^ — a*a) — D{a), the displacement operator, and £ is 
that of the RWA, then 

C'{t) = C{t) + (^{lil + X)<j + (at , _ , at) _ + X)a* + (a . - . a). (29) 

It is easy to see that if 

cr(i) = cr(0)exp(-ir2(t) - A(i)), cr*(i) = (cr(t))*, (30) 

then both p and p' satisfy the same equation. That is, we have shown that D{a{t))p{t)D^ {a{t)) satisfies Eq. 
( pT| ) whenever p{t) does the same. We remark that this result does not depend on the temperature of the 
bath. 

We now turn to the evaluation of the density matrix evolved with the superoperator found above. We 
chose initial states relevant from the point of view of quantum optics. As a first initial state we choose the 
system's ground state, p(0) =| 0)(0 |. Since 

e^^ I 0)(0 I = I 0)(0 |= e'^ \ 0)(0 H 0)(0 |, and (31) 

oo ^ oo 

e-^|0)(0| = ^^(at)"|0)(0|a" = ^:."|n)(n|, (32) 



we have 



oo ^ oo 

Pit) = U{t) I 0)(0 h E YTN + ^'^^^ ' 1= E^"W I I ■ (33) 



The above formula displays the so called decomposition in natural orbits where the quantities P„ can 
be directly interpreted as probabilities. We can write the evolved density matrix in the alternative form 
p{t) — exp((l + l/M)a^a)/{\ + M), which is the form of an instantaneous thermal density matrix, with 
Mit) =< a^a > (t). Had we chosen an initial thermal state, with mean number of excitations n(0), the 
density matrix would have remained a thermal state, but now M{t) = n{0) exp(— 2A(i)) +A/'(t). If we use the 
instantaneous oscillator frequency = a; + 5, it is possible to define an instantaneous temperature through 
the relation T{t) = hiuj + d) / [Ub lii(l + l/A/')), with ks the Boltzmann constant. Moreover, we have obtained 
a physical interpretation for the quantity M{t): it is the mean number of excitations of the main oscillator 
at time t when it was initially prepared in its ground state. 
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To calculate the density matrix for an initial Fock state it is better to write the evolution superoperator 
in the form U{t) = v exp(wTZ) exp{z' J^) exp(xA4) exp(yV), where w,x,y are given in (p^), and z'{t) — 
(l+A/')exp(2A). We use 

m I 

e^-^e^-^e^^ I m)(TO I = ^(e^+y)"'-fc(^e^+y)fc \jn~k){m-k\ 

m I 

= E7;;^^W(e^^')'(-e^^')™"'|fc)(^l' (34) 

k=0 ^ ' 

e-^\m){m\ = T^j^^lkHkl (35) 
to see that the density matrix at time t is given by p{t) = X]fc°=o ^ni,s(^) I (s li with 

»! (m -*)!(»-(■)! I, ■ ' ' 

Since the former density matrix have been expressed in terms of natural orbits, the quantities Pm,s{t) are 
readily interpreted as probabilities. The transformation property discussed above allows us to write the 
evolution of an initial generalized coherent state | am) = D{a) \ m), where D is the displacement operator 
and I n) the n-th number state. We have U{t) \ aQm){(To'm |= '^'^^q Pm,s{t) \ (j{t)s){a{t)s \, with Pm,s{t) 
givenby (H) anda(i) by (|^. 

One interesting point to be investigated is if there exists an asymptotic density operator. Provided 
that our environment is such that limijoo A(<) oo and limtjoo A/'(t) — rioo, the asymptotic evolution 
superoperator can be written as 

1 77 

limW(t) ^ — exp(— 2^^7^)(| 0)(0 | •)exp(J)(. | 0)(0 |), (37) 

tjoo i + rioo J- + 



which applied to a generic normalized initial density p(0) gives 

1 TL 

Poo = — exp(^)p(0)(0|exp(— 2^7^) |0)(0 

l+Uoo l+rioo 



^ (Trp(0))exp[(l + — )ata]|0)(Oh-^ exp[(l + —)«+«] | Q) (0 | . (38) 



1 -\- Tloo ^oo 1 ^oo 

Thus, whenever the established conditions are met, the density operator approaches asymptotically to a 
thermal state with a mean number of excitations equal to that of the environment. The existence of a unique 
asymptotic density can not be taken for granted: in the model of decoherence without damping studied in 
references even when the coefficient of decoherence grows indefinitely with time, the asymptotic state 
depends on the initial state. 

The normal order characteristic functionalC*^"-' (^, t) given by 

C(")(^,$t^i) ^ Tre'^^^^^'XO = Tre*«'^'e'«*"Z^(Op(0) = TrZY'f (t)e'«"^*«*V(0), (39) 

where t{{t) is the evolution superoperator and W{t) its adjoint, is the generating functional of the normally 
ordered moments. After a somewhat lengthy but straightforward calculation we obtain 

— g-2AA5a;itrp^gi^oxp(-A+in)at^i5* exp(-A-in)a^^Q-j^ j-^q-j 

If we calculate (a) (t) and (a f a) (t) using 

m) - k=o=«t, {aU)it) = a^C'(")(^,e^t) lc=o=?t, (41) 

we arrive at the same results as before. These can be compared to those obtained in reference |q]. 
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4 The Zero Temperature Limit 

The zero temperature limit has its own special interest, both as an approximation at low temperatures, and 
as the relevant case for leaking Bose-Einstein condensates. In this case the evolution superoperator can be 
expressed as U{t) = e'^^e^^'e^--'^ with 

S:{t) = -in{t)-A{t) = y*{t), = 1 - exp(-2A(t)). (42) 

Since the vacuum is solution of the master equation, the transformation property discussed above indicates 
that 

Dia{tm{0\DH<^it)) = W{t)){a{t)\ ^ U it)\ao) {ao\ (43) 

also solves the master equation. That means that initial coherent states evolve preserving their coherence, 
not matter what the details of the environment. A measurement of the norm and the phase of an initial 
coherent state is enough to determine the functions Q and A, and hence, in principle, to determine the 
evolution of any other initial state. 

The evolution of an initial Fock state is also easily calculated. 

m I 

mH{m\^J2p,,Utm{k\^p,,^{t) = fc), (e-^-^'^^)^(l - e-^^W)"-''|fc)(fc|. (44) 

From this solution we can generate another solution: if we have an initial generalized coherent state, |m; uo) = 
D(ao)\'m), it evolves into a mixture of generalized coherent states, as given by 

m 

Uit)\m;ao){m;aa\ = "^pk,Ut)\k;<j{t)){k;a{t)\. (45) 

k=0 

Applying the evolution operator Z//(i) to an initial density matrix element |cro)(cr'o|, one obtains 

mWo){<T'o\ - /%ttM t))i^'it)\- (46) 
{a'{t)\cr{t)) 

Notice that an heuristic argumentation also leads to (|4^). Indeed, since a coherent state remains coherent, 
we expect U {t)\<To) {(j' o\ = N{t)\<T{t)){a' {t)\. As long as the exact dynamical equation for p preserves the 
trace, ■§;Ttps{t) = Tr(^^^^|^) = TT{C{t)ps{t)) — 0, the normalization factor N{t) cannot be other than that 
of Eq. (^6|). The evolution of an initial even [paoe) or odd cat [paoo) state can be calculated from ( ^6| ) and 
(H). Indeed, 

P.oe(o) = N,(o){oo){\oo), I - ao)) ) ) , (47) 

where Afe(o)(o'o) = (1 + {~){^'^oWo))^^ /"^ is a normalization factor, evolves as follows 



(48) 



(1 (-)(-To|ffo) \ / I \ \ 
iMSS'\ I a.) • 

We can rewrite Eq. (|4^) in a more convenient way, in terms of natural orbitals, as 

{t)Paoe(o) = Pe(o) (OP<T(t)e(o) + (OPrT(t)o(e) , (49) 

with 



e(o),. ^ l l + {-){-a{t)\a{t)) ( (-aolao 

2 l + (-)(-aoko) ^ {~a{t)\a{i)) 



(50) 



^^(°)^'^ - 2 l + (-)(-a„|ao) V {-<r{t)\a{t)) ) ■ ^^^^ 
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Observe that an initial cat state evolves as a mixture of even and odd cat states. Eq. ( |50| ) gives the 
probability of the cat state of the same parity as the initial state, and Eq. the probability of the cat 
state of the other parity. 

The evolution of an initial squeezed state is also easily computed. We notice that the density operator 
for the squeezed vacuum /3(C), can be written as 

p{C) = lim p(C,7) = lim (1 - e-^)S{Oe-''''''' S\C), (52) 

7 — ^00 7 — >oo 

where 5(C) — exp((C(at)^ — ("l'a^)/4). Setting ( = ^expicf), the following expressions for the second moments 
of p{C, 7) are found 

((at)2) = ^sinh(e), (a^) = f^sinh(e), ({a, at}) = cosh(0. (53) 
The system of equations for the second moments, 

y.-2i(.n + A,(0, (54) 
^^=2^(,n-A)((at)2), (55) 

^i^ = -2-({a,at}) + _, (56) 



is integrated to give 



((a2)(t))(t) =cxp(-(-)2zl7-2A)((a2)(t))(0), (57) 
{{a,a)}){t) = exp(-2A)({a,at})(0) + (1 - exp(-2A))/2. (58) 



From these relationships we find U{t)p{() — p{({t),'j{t)), with 

at) = at) exp{i(j)it)) (j}{t) =<j)o-2n (59) 

at) = ArcTanh( 7) (60) 

' ^cosh(^o) +exp(2A) - 1^ ^ ' 

7(<) = ArcCoth(Y/e-4A + 2 cosh(^o)e-2A(i _ g-^A) + (1 _ q-^^)^). (61) 

Now the evolution of a general initial state p((T, Q — D{a)p{C,)D'^ [a) is 

U{t)p{aX) = piait),at),lit))=Diait))piat),lit))D^{ait)), or (62) 
Zi(t)(|ao,Co,n = 0)(ao,Co,n = 0|) = i?(a(t))5(C(i))(l - e-^W)e-^W'^'^5t(C(t))i^t (^(i))^ (53) 



where a-(t),C(i) and ^(t) are those of Eqs. and (p9[-|6lD, and 

\aoXo,n)=D{<7{t))S{at))\n). (64) 

The corresponding expansion in terms of natural orbits is 

00 

w(i)(ko,Co,o)(tTo,Co,o|) = ^p„(<)k(0,C(0,^)('^(i),C(i),'^l 

00 

= ^e-"^(*)(l-e-^W)|a(t),C(t),n)(a(t),C(t),n|. (65) 



n=0 
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5 Determination of the Master Equation Coefficients 



Let's return to Eq.(|10[). If we call the exact eigenfrequencies of (|1|), we can write (remenber that greek 
indices can assume the value 0, while latin indices do not) 

a,it) = J2 U;,U^^e-'^-*a„{0) = ^ Z,,a<,(0). (66) 

For the sake of convenience we write 

'7 = ^00, Ik = Zok, = Zko, Tki = Zki- (67) 

Using (|6|) in Eqs. (0) and (p^), we obtain the following expressions for rj^ and 7^;, 

Afe Ak-fi , , 

Vk = — Iki = Tfci . (68) 

V V 

From (|l|) and ^ we obtain 

A = - ^ CkMvk), A' = - ^ Cfe(2nfe(/3) + l)Im(/3fe) = - ^ Cfe(2nfe(/3) + l)Im(^ 7,7^,). (69) 

Since J7 is unitary we have 

E^'^^^ = E^;oC^MAC/.fet^:Ae"^("'''-"'^)*"E^;ot^MoC/.feC/;oe-(-""""')* 

= J2KoU.ke-'^''^-'^^'S^,-rjAl^Y.U;,U,k-V^l=-V^*k, (70) 

and 

E ^'^r = E C/;oC/^AC/.oC/;;,e-*("'--'^)* - E C/;oC/^ot/.oC/:oe-'("'^^"'-^* = 1 - W*- (71) 

I \fJ,U fiv 

Using (^ and (|7l|) we get 

= E ^'T^*^ = E ^'^H - ^ E ^'7' - ^^fc - - W*) - -ril- (72) 

Thus we finally notice that \'{t) can be expressed as a sum 

A' (t) --Ecfc(2nfc(/3) + l)Im(%) =A(t)- 2 EcfcrH.(/3)Im(7?fe) -A(t) + £(<;/?). (73) 
fc fe 

Observe that lim^_»oo e(t; /3) = 0. This proves that indeed we have the equality \{t) = X'{t) at zero temeper- 
ature, as can be seen by looking at their definitions as given in Eq. (|69|). 

We have seen that S, A and e (or il, A and Af) are all we need to characterize the effect of the bath on 
the main oscillator. The Heisenberg equations of motion, 



d = -iwoa - iE^fca/c, (74) 

k 

dk = -iuJkak - icfca, (75) 
can be solved in a number of ways. For example, an implicit method gives 

a(i) = a(0)e~*"(*)-'^(*) -zEcfcafc(O) / dTe'"^''^'-'''>e-'^'^^^~^'^^\ (76) 
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where r](t) satisfies the integrodiff'erential equation 



f] + iuorj + 



f drY,cle^-^^'--^r,{r)^0, (77) 
Jo 



subject to the initial condition etap ~ 1. Using the %(t) impUcitly defined above in eq (|T^, and taking into 
account the equation satidfied by 77, we obtain \ + iS = —iuj + d{\nri)/dt, or ri{t) — exp(— A(t) — iil{t). We 
also find 

e{t) - ^^j- e^^W E I / dre-^-^(*--)e-*"(^)-^(^) P nfe(/3) . (78) 



dt , 

\ k 







Had we used normal modes wewould have arrived at the expressions 



= E ^1 ' and (79) 

e{t) = Re — > — ^iZ-i^ . 80) 

We moreover would have noticed that, in order not to have inverted oscillators the following condition has 
to be fulfilled, cv > cl/uk- Since ( [tt] ) can be hard to solve, methods to obtain approximate solutions are 
welcome. For instance, we can expand in powers of Ck to second order, a(t) — exp(fco + cA;i + c^fc2), to obtain 

- in{t) - A(i) = -iujot -J2^k [ dti [ dtae"'^"'"^"'')*^ (81) 

k "^0 

We expect (|8^) to be valid for small Ck- It can be shown that this is the time dependent Born-Markov 
approximation, which is valid also for "strong" coupling and short times [l^. 



6 Experimental Characterization 

We observe that the calculation of the mean energy, and of the entropy of the above examples of initial 
conditions only require the knowledge of A(t). However, if for example one wants to measure the Wigner 
function of any field density matrix, one would need to know n{t) as well. To determine this function 
one takes advantage of the experimental setup to measure the Wigner function^, in which a two level 
atom is prepared in the excited state |e), sent trough an array of two low-Q cavities Rl and R2 and one 
high-Q cavity C between them, and is detected eventually. The field in C is displaced by the operator 
D{a) = exp(aa^ — a*a), by a microwave source connected to it. Rl and R2 behave as "rotation" operators 
in the Hilbert space of atomic states, | e) ^ (| e) + | g))/V2, and | g) (-6"*^ | e)+ | g))/V2, with 
^ = in Rl and ^ = 7r/2 in R2. The dispersive atom-field interaction in C produces entanglement: The 
field component associated with | e) suffers the action of the operator exp{iTr{a^a + l)/2), while the one 
associated with | e) suffers the action of exp(— i7ra^a/2). It was shown in Ref. Q that for this experimental 
arrangement 

AP = P,-Pg = W{^a,-a*,t)/2, (82) 

where Pe {Pg) is the probability to detect the probe atom in the upper (lower) state |e) {\g)), and W{—a, —a* , t) 
is the value of the Wigner function corresponding to the field in thehigh-Q cavity at the time t the probe 
atom exits this cavity. Notice that, for the sake of convenience the normalization of W has been changed by 
a factor of tt. The Wigner function of the coherent state ( ^ ) is 

Wia,a*,t) = 2e-'("(*)-")("(*)-")'. (83) 
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Let's notice that W{a, a* ,t) presents a maximum at a = a{t) — (Toexp(— A(i) — iCl{t)). Since A can 
be determined from a photocounting experiment, we just need to measure Q{t). We choose a = a(t) = 
ctq exp(— A(<)) exp(— and adjust to maximize AP. If = f7(i)mod(27r), then we obtain 

AP = 1, otherwise, it wiU be smaller than one, AP = exp(— Suq exp(— 2A) sin^((ri — <i>)/2)). Choosing 
different exit times t we can map fl{t). Notice that what we do, in fact, is find the right a{t) that brings 
the field at the cavity C to the vacuum. If we succeed in so doing the probe atom is rotated in Rl, getting 
in a superposition of upper and lower states, flies free till R2 where the rotation is reversed, and goes to the 
detectors again in the upper state. 

7 Concluding Remarks 

In conclusion, we have shown that, if the optical cavities can be modelled through the hamiltonian (|l|), then 
only three, experimentally measurable real functions are necessary to characterize their behavior: the mean 
photon number of an initial ground state, the instantaneous frequency and the rate of change of the number 
of excitations. Measuring thes quantities once will enable one to get the Wigner function for any initial 
condition. So, provided the adequacy of the RWA model has been tested the present work provides for an 
alternative way to construct the time evolution of Wigner's functions. 
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